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Résumé : 

L’identification de paramètres d’une loi élastique et de deux lois de comportement élasto- 
plastique à écrouissage cinématique intégrée à la corrélation d’images (CINI) est menée sur une 
éprouvette biaxiale cruciforme. La géométrie est optimisée spécifiquement pour l’identification 
d’une loi élastique linéaire. Cette optimisation s’exprime dans le formalisme de la mesure des 
champs de déplacement et considère l’ensemble de l’essai depuis l’acquisition d’images et la me¬ 
sure d’efforts jusqu’à l’estimation des paramètres constitutifs. Par conséquent, l’ensemble des 
incertitudes existantes est pris en compte dans la définition du protocole expérimental. 

Abstract : 

Integrated digital image corrélation (IDIC) is applied to identify the constitutive parameters 
of an elastic and two elastoplastic laws with kinematic hardening. An experiment is conducted 
on a cruciform specimen in a biaxial setup. The geometry of the specimen is optimized to 
allow for the least uncertainty in the identification of parameters of linear elastic law. The 
optimization takes into account the complété experimental and identification processes. Thus, 
full held measurements, direct numerical simulations and identification procedures accounting 
for uncertainties, ail contribute at their best to the resulting experimental protocol. 

Mots clefs : Corrélation d’images, Elastoplasticité, Identification 

1 Introduction 

L’identification rapide et robuste des propriétés mécaniques sous chargements monotone et cy¬ 
clique est un enjeu majeur pour la conception mécanique. Par exemple, l’objectif de réaliser 
des structures plus légères et plus sûres, conçues dans des délais toujours plus courts, nécessite 
une description plus précise et plus large des comportements des matériaux retenus. Un des 
premiers exemples concerne la caractérisation du comportement mécanique de tôles d’acier dont 
l’orientation de la structure vis-à-vis de la direction de laminage [T] affecte le comportement. 
Les enjeux sont d’identifier les lois et les paramètres de manière robuste et rapide. La technique 
d’identification intégrée à la corrélation d’images (CINI) apparaît être une solution capable de 
répondre à ces objectifs. Cette dernière a déjà été mise en œuvre pour la caractérisation de 
différents comportements mécaniques dllIlIS] et fera l’objet d’une brève présentation. Néan¬ 
moins et afin de tirer le meilleur parti de celle-ci, une technique d’optimisation est proposée et 
dont une application est faite pour les paramètres d’un loi élastique linéaire et isotrope sur une 
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éprouvette cruciforme. Enfin, l’identification des paramètres d’une loi linéaire élastique et de 
deux lois élastoplastiques à écrouissage cinématique, i) linéaire et ii) non-linéaire est réalisée. 

2 Outils 


La corrélation d’images est conduite entre une image de la configuration de référence carac¬ 
térisée par son niveau de gris sur chaque pixel f{x) et une série d’images dans des configurations 
déformées g{x) [ül[7]. L’hypothèse de la conservation des niveaux de gris au cours du temps per¬ 
met de poser comme seules inconnues du problème les composantes du champ de déplacement 
{u{x)) entre ces images. Le problème à résoudre s’exprime par la minimisation d’une fonction¬ 
nelle Xj = \/Nt X/(^) dui est la différence quadratique moyenne sur l’ensemble de la région 
étudiée entre les images f et g corrigées du champ de déplacement mesuré u{x,t) 


Xf{t) = 2^jNQ + u(x, t), t) - f{x)f 


( 1 ) 


où Nçi est le nombre de pixels dans la région sur laquelle la corrélation d’images est menée, et 
Nt le nombre d’images acquises au cours de l’essai. Pour l’identification intégrée, les inconnues 
sont les paramètres inconnus de la loi de comportement choisie. Le problème revient à trouver 
le meilleur jeu de paramètres qui minimise globalement la fonctionnelle x^f rapport aux 
champs de déplacement appartenant au sous-espace engendré par les champs de sensibilité aux 
paramètres 


{ 5 p} = 


( 2 ) 


où [M]jr)ic = 1 /{2xJ)[SuY'[M][Su\ est le hessien cinématique, [M] la matrice de corrélation 
d’images pour une discrétisation choisie, [S'cr] est la matrice contenant toutes les composantes 
des champs de sensibilité [5] vis-à-vis des paramètres recherchés ([5'j/] = \du/dp]), {p} est le 
vecteur contenant les paramètres recherchés et enfin {6} = (/ — g){^f). Par ailleurs, les efforts 
mesurés peuvent aussi servir à l’identification. L’écart entre les valeurs mesurées et simulées 
s’exprime dans la fonctionnelle suivante que l’on cherche à minimiser 


Xf = 


1 


-{Fm - F,Y{F^ - FJ 


( 3 ) 


NpNtil 

où N F est le nombre de données mesurées sur les efforts à chaque pas de temps, et Nt le 
nombre de temps de mesure. Pour cette seule fonctionnelle, la linéarisation du problème revient 
à résoudre successivement jusqu’à convergence 

1 


{ôp} = ^[H]Ÿ[SFf{Fm-F,} 


( 4 ) 


où \H]f = I/xfISeYISf] est le hessien statique, [S'f] est la matrice contenant les sensibilités en 
effort, {Fm} et {FJ sont respectivement les efforts mesurés et calculés aux différents instants 
considérés. 

Les deux types de mesure sont combinés en normant chacune des fonctionnelles par l’ampli¬ 
tude de leurs incertitudes respectives. Le poids spécifique provient d’une formulation bayésienne. 
L’hypothèse sous-jacente est l’existence d’un bruit blanc gaussien, i.e., non corrélé dans le temps. 
Ainsi, si seul le bruit est présent (aucune erreur de modèle) la valeur attendue de x^ est 1. Cette 
dernière s’exprime par 

2 _ Nq 2 , 2 


X = 


+ Nn + Np' 

La matrice de covariance des paramètres identifiés s’écrit alors {[H] = [M]jfijc + [H]f) 


-Xf 


( 5 ) 




( 6 ) 
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3 Loi de comportement 

L’identification est réalisée pour une loi élastique linéaire (^) et deux lois élastoplastiques, 
la première {B) à écrouissage cinématique linéaire et la seconde (C) à écrouissage cinématique 
de forme exponentielle. Plusieurs auteurs ont déjà proposé des études sur l’identification des 
paramètres de loi de comportement élastoplastique 00 . La déformation totale est séparée en 
une contribution élastique et plastique etot = CeZ + ^pi où la loi élastique est isotrope et linéaire 
telle que Cei = i^cr — ■gtr(cr)J. v et E sont respectivement le coefficient de Poisson et le module 
d’Young et I est le tenseur identité. Le seuil de plasticité est défini par la contrainte équivalente 
de Von Mises 

/(o- - X) = - Xy : (cr - X)'') - ctq = 0 (7) 

où (To est la limite d’élasticité macroscopique du matériau et — l/3tr(*)/ est la partie 

déviatorique du tenseur •. Dans le cadre d’une modélisation linéaire de l’écrouissage, la variable 
d’écrouissage cinématique X s’exprime par X = ^Cépi. Pour la seconde loi à écrouissage ci¬ 
nématique non-linéaire de forme exponentielle, son expression est X = ^Cépi — cXp où p est 
le taux de déformation plastique cumulée [101 111] . Enfin, les paramètres sont exprimés suivant 
une échelle logarithmique qui permet de privilégier les incertitudes relatives. La paramétrisa- 
tion retenue s’exprime selon {q] = (log({p/po}) où {p} est le jeu des paramètres recherchés 

M = 

4 Optimisation 

Il n’est pas aisé de formuler un critère d’optimisation universel lorsque plusieurs paramètres 
sont en jeu. Dans le cas ici présenté, l’optimisation est naturellement multiparamétrique et un 
critère doit être formulé. La stratégie choisie est de viser l’incertitude la plus faible pour le 
paramètre le plus sensible au bruit. L’optimisation consiste donc à minimiser la plus grande des 
valeurs propres associées à [Cp\ afin de minimiser l’influence du bruit sur les paramètres associés. 
Ceci est équivalent à maximiser la plus faible valeur propre de [H]. Les paramètres matériaux 
correspondent à un acier inoxydable à durcissement structural 17-7 PH [T^] {E = 200 GPa, 
V = 0.3 et (To = 1300 MPa). A cette étape, une hypothèse de champ moyen est exploitée du fait 
de la méconnaissance a priori du mouchetis. Elle permet d’exprimer la matrice de corrélation 
comme 

ipi{x) ■ ipj {x) (8) 

n 

où Gj = ((V/)^) est le gradient de niveau de gris moyen et sont les fonctions de 

forme choisies. L’incertitude des efforts provient de deux sources indépendantes, la première 
est l’incertitude des capteurs de force (dépendant du niveau d’effort). La seconde concerne les 
conditions aux limites de type Dirichlet (obtenues à partir d’une solution initiale du champ 
de déplacement à l’aide de la CIN). En effet, des niveaux de déplacement affectés du bruit 
d’acquisition vont conduire à des efforts simulés erronés. Le modèle d’incertitude global retenu 
s’exprime selon 

(9) 

où Pi est l’écart type de l’incertitude qui dépend de l’effort imposé et po une valeur constante 
correspondant à l’erreur due aux conditions aux limites imposées issues de la CIN. Cette der¬ 
nière est évaluée à chaque itération et tient compte de l’évolution des valeurs des paramètres. 
Le tableau [T] réunit les valeurs d’incertitudes retenues. L’historique de chargement est non pro¬ 
portionnel et impose dans un premier temps une amplitude de déplacement di = d* suivant 
la direction ei en maintenant le déplacement suivant 62 à ^2 = 0. Une fois celle-ci atteinte, le 
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déplacement est imposé suivant jusqu’à la même amplitude d 2 = d*. Puis la décharge est 
équibiaxiale jusqu’à di = d 2 = 0. Dans le plan (di, ^ 2)1 Is trajet suit la forme d’un triangle, d’où 
sa dénomination « chargement triangle ». 


Table 1 - Résolution et incertitude des déplacements et efforts mesurés. Les valeurs sont éva¬ 
luées à partir de cas expérimentaux. La taille du pixel correspond à a = 13.5 pm 


Quantité 

PO 

Pi 

If 

Gf 

a 

Valeur 

2.5 N 

2.9x10"^ 

280 niveaux de gris 

3800 niveaux de gris/pixel 

13.5 pm 


Le paramètre géométrique optimisé est le rayon du congé de raccordement reliant les bras de 
l’éprouvette. La hgure[ÏJa) présente une des géométries ainsi que le maillage associé. Le rayon 
de raccordement présenté Figure [I^a) est r = 2 mm. Figure mb) montre les valeurs propres 
de la matrice hessienne [H] en fonction du rayon de raccordement évalué pour l’ensemble de 
l’historique de chargement. L’amplitude du déplacement imposé est adaptée à chaque valeur 
de r de telle sorte que la valeur maximale de la contrainte équivalente de Von Mises soit égale 
à la limite d’élasticité macroscopique du matériau. La valeur optimale du rayon est évaluée à 
r = 1.7 mm, correspondant à la valeur maximale de A /_2 (figure [ÏJb)). 



Dimension (mm) 



Figure 1 - (a) Maillage et géométrie de l’éprouvette biaxiale à optimiser pour un rayon 
de raccordement égal à r = 2 mm. (b) Valeur propre A/ de \H] pour le jeu de paramètre 

(iog(F;/F;o),iog(./.o)) 


5 Identification 

L’identification est réalisée sur un essai biaxial à efforts imposés et à plusieurs cycles à l’aide 
du dispositif expérimental présenté Figure [21 Le matériau est l’acier inoxydable 17-7 PH avant 
précipitation par traitement thermique (le comportement est plus ductile). Ces caractéristiques 
sont différentes de celles retenues pour l’optimisation avec une déformation totale potentielle plus 
importante (e « 35%) mais avec une limite élastique plus faible (ag = 300 MPa). Le chargement 
est aussi non proportionnel et triangulaire mais est à efforts imposés. Les images sont acquises 
à l’aide d’une caméra 16-bits PCO Edge et un objectif télécentrique de grossissement xO.5. 
Les paramètres initiaux sont ceux identifiés dans un premier temps uniquement à partir des 
efforts de réaction (Equation [3]) . La Figure llfb) montre les résidus de corrélation à convergence 
issus du calcul CIN initial (noir) et des calculs CINI (couleurs). Après un début d’essai pour 
lequel les résidus CIN et CINI sont sensiblement identiques, ceux-ci s’écartent une fois que les 
déformations deviennent plus importantes. Le résidu pour la loi A s’écarte de la solution CIN 
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Figure 2 - (a) machine de traction compression biaxiale et (b) résidus de corrélation CIN (noir) 
et CINI (color). 


à partir du pas de temps t = 25, tandis que les résidus de i? et C s’écartent à partir du pas 
de temps t = 80. Il apparaît à la fin de l’expérience que le résidu d’identification le plus faible 
est celui obtenu pour la loi C. Enfin, une fonctionnelle de régularisation {xr) est associée à la 
fonctionnelle à minimizer (x) afin de prévenir une évolution irraisonné des paramètres dont la 
sensibilité est inférieure au bruit. XR est une fonction convexe et son minimum est égal zéro. 
Celle-ci s’exprime 

xl = {q}[C^V{q} ( 10 ) 

où [C^] est la matrice de covariance associée à la résolution limite choisie. Enfin et puisque XR 
et X sont sans dimensions et dimensionnés à 1 la fonction à minimiser régulariser devient 

([Ef]+A*[/]){<5g} = {b} + A*M (11) 

où A* est le paramètre de régularisation. Celui-ci est choisi de tel sorte que les paramètres dont 
l’influence sur les grandeurs mesurables soient inférieurs à l’incertitude de mesure ne subissent 
pas d’évolution. La table 2 montre les paramètres identifiés sur le premier cycle où un très 
faible niveau de déformation plastique est atteinte et où celui-ci n’est pas suffisant pour révéler 
l’influence des paramètres matériaux. Par conséquence et grâce à la régularisation leurs valeurs 
ne subissent pas de variation. 


Table 2 - Identified parameters and identification residuals for the three laws for the first cycle 
( 7 p is expressed in % oî p). 


loi 

X 

X/ 

XF 

E (GPa) 
200 

JE 

U 

0.3 

'ï<' 

cro(MPa) 

300 

7o-o 

CfGPa) 

10 

7c C 7c 

10 

A 

2.10 

2.10 

15 

174 

0.4 

0.31 

0.09 

— 

— 

— 

— — — 

B 

2.09 

2.09 

8.9 

195 

0.21 

0.297 

0.1 

306 

0.25 

10 

— — — 

C 

2.09 

2.09 

8.9 

195 

0.25 

0.297 

0.1 

306 

0.3 

10 

— 10 — 


La Figure [3] montre les évolutions de ceux-ci au fur et à mesure des itérations. Pour la loi 
A, le paramètre ly converge vers v = 0.50, en effet le comportement étant plastique, il apparaît 
quasi-incompressible. Les autres paramètres convergent relativement rapidement pour les lois 
B et C. Le tableau 3 montre les résultats de l’identification sur l’essai complet. Le module 
d’Young identifié pour la loi A n’est pas physique. En effet, la déformation étant principalement 
plastique, et la loi élastique linéaire, le module d’Young identifié est ajusté à une valeur trop basse 
pour rester compatible avec les efforts mesurés. Une loi élastoplastique est plus satisfaisante, et 
conduit à un résidu d’identification plus faible. Pour la loi à écrouissage cinématique non-linéaire. 
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-E/Eç, 

-vjva 

-CICo 

- o-/cro 

-7/70 


(a) 


(b) 


(c) 


Figure 3 - Evolution des paramètres matériaux respectivement pour les lois (a) A , (b) B et 
(c) C. 


ce résidu est de l’ordre de 6.5 fois le niveau de bruit, montrant encore une erreur de modèle, 
mais cependant une amélioration sensible par rapport au cas élastique. Le résidu est légèrement 
plus élevé pour la loi à écrouissage cinématique linéaire. 


Table 3 - Identified parameters and identification residuals for the three laws for the whole 


loading history 


loi 

X 

Xf 

XF 

E (GPa) 
200 

lE 

V 

0.3 

Xv 

f7o(MPa) 

300 

7cro 

C(GPa) 

10 

le 

C 

10 

7c 

A 

15 

14.9 

1820 

8.8 

0.25 

0.499 

0.002 

— 

— 

— 

— 

—■ 

-- 

B 

7.1 

7.1 

113 

157 

0.15 

0.31 

0.04 

480 

0.15 

6.7 

0.15 

— 

— 

C 

6.3 

6.3 

100 

148 

0.16 

0.30 

0.04 

423 

0.16 

8.6 

0.16 

10.8 

0.03 



(a) Loi A 


(b) Loi B 


(c) Loi C 


(d) GIN 


Figure 4 - (a), (b) and (c) différence entre les champs de déplacements issues de CIN et de CINI 
pour les lois T, B et C à la fin de l’expérience, (d) composante Ui du champ de déplacement 
CIN au même instant 


La figure |4] présente la composante Ui du résidu champ de déplacement à la fin de l’expé¬ 
rience, i.e., t = 350 s, où encore Cmax = 18.9%. De la gauche vers la droite, le résidu pour la loi 
A, i.e., élastique linéaire, la loi B et enfin C. Pour la première, le résidu correspond directement 
aux déformations permanentes qui ne sont pas prédits par une loi linéaire élastique. Concernant 
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les deux autres composantes du déplacement résiduel, la loi B est à un niveau faible (dynamique 
de déplacement réduite d’un facteur 2) et la loi C améliore encore ce niveau de résidus. Malgré 
tout, des déformations (hors de l’incertitude de mesure) restent inexpliquées. Il y a donc ma¬ 
tière à généraliser la forme algébrique de la loi de comportement pour mieux s’approcher des 
observations expérimentales. 

6 Conclusion 

Une nouvelle approche d’identification optimisée a été proposée. Sa mise en œuvre sur une 
éprouvette biaxiale a permis d’identifier une géométrie minimisant l’influence des incertitudes 
de mesure sur les paramètres à identifier dans le cas d’une loi élastique linéaire. Pour cela, 
l’ensemble de la chaîne de mesures exploitée pour l’identification intégrée est considéré. L’op¬ 
timisation est spécifiquement réalisée pour la loi dont on cherche à identifier les paramètres et 
dont on suppose a priori le comportement expérimental. L’étude expérimentale est cependant 
effectuée dans le cas d’un trajet de chargement qui explore largement le domaine plastique. 
L’identification a été réalisée pour trois lois de comportement et les paramètres ont été identifiés 
avec un niveau d’incertitude faible. Néanmoins, le résidu d’identification, relativement élevé en 
comparaison du résidu GIN, montre clairement l’existence d’une erreur de modèle résiduelle, très 
forte pour une loi élastique linéaire, et bien moindre pour une loi élasto-plastique. Ceci permet 
de classer les modèles en fonction de leurs qualités respectives vis-à-vis des données expérimen¬ 
tales. Dans la présente étude quatre indicateurs d’erreur ont été considérés, 

- Quand les analyses GIN et CINI sont employées, l’écart par rapport à l’hypothèse de 
conservation de niveau de gris est évalué. Il estime la qualité de mesure lorsque les images 
déformées sont corrigées par le champ de déplacement mesuré et comparées à l’image de 
référence. 

- GIN et CINI peuvent aussi être comparées en calculant les déplacements mesurés par les 
deux approches. Lorsque le même maillage est employé, la comparaison est directe. 

- Le résidu d’effort et un autre moyen de caractériser la qualité de l’identification. Ici, FEMU- 
F et CINI peuvent être comparées car les deux approches inclues les efforts dans leurs 
formulations respectives. 

- Finalement, la méthode d’identification intégrée CINI possède un indicateur de qualité qui 
considère toutes les sources d’informations utilisées dans un unique formalisme grâce à la 
normalisation par les variances des l’incertitudes correspondantes. 


L’ensemble de ces indicateurs permet à l’utilisateur d’aboutir à une identification globale 
mais aussi avec évaluer la qualité de l’identification et l’erreur de modèle (i.e., choix d’une loi de 
comportement, discrétisation en éléments finis). Cette étude cherche à présenter et à construire 
une route pour la compréhension des phénomènes mécaniques, la conception optimisée des essais 
expérimentaux et l’identification pour valider mais aussi invalider les modèles étudiés. Dans des 
études futures, il est envisagé d’étendre cette technique à des matériaux différents tels que 
des polymères (comportement hyper-élastique), céramiques (comportement fragile) ou encore 
composites (comportement anisotrope). Enfin, des problèmes à l’échelle microscopique avec la 
prise en compte de la microstructure pourront être envisagés. 
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